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Abstract. The analysis of the broad iron line profile in the X-ray spectra of 
active galactic nuclei and black hole X-ray binaries allows us to constrain the spin 
parameter of the black hole. We compare the constraints on the spin value for two 
X-ray sources, MCG-6-30-15 and GX 339-4, with a broad iron line using present 
relativistic line models in XSPEC — LAOR and kyrline. The LAOR model has 
the spin value set to the extremal value a = 0.9982, while the kyrline model 
enables direct fitting of the spin parameter. The spin value is constrained mainly 
by the lower boundary of the broad line, which depends on the inner boundary 
of the disc emission where the gravitational redshift is maximal. The position 
of the inner disc boundary is usually identified with the marginally stable orbit 
which is related to the spin value. In this way the LAOR model can be used to 
estimate the spin value. We investigate the consistency of the LAOR and kyrline 
models. We find that the spin values evaluated by both models agree within the 
general uncertainties when applied on the current data. However, the results are 
apparently distinguishable for higher quality data, such as those simulated for the 
International X-ray Observatory (IXO) mission. We find that the LAOR model 
tends to overestimate the spin value and furthermore, it has insufficient resolu- 
tion which affects the correct determination of the high-energy edge of the broad line. 



Introduction 

The broad emission iron lines are well-known features found in about two dozens of spectra of active 
galactic nuclei and black hole binaries. They are supposed to originate close to the black hole by the 
reflection of the primary radiation on the accretion disc. The spin of the black hole plays an important 
role in the forming of the line shape. Especially, it determines the position of the marginally stable orbit 
which is supposed to confine the inner edge of the accretion disc (see Figure [TJ . The innermost stable 
orbit occurs closer to a black hole with a higher spin value. However, the spin affects also the overall 
shape of the line. 

Over almost two decades the most widely used model of the relativistic disc spectral line has been 
the one by Laor [1991] , which includes the effects of a maximally rotating Kerr black hole. In other words, 
the LAOR model sets the dimensionless angular momentum a to the canonical value of a = 0.9982 - so 
that it cannot be subject of the data fitting procedure. Dovciak et al. [2004] have relaxed this limitation 
and allowed a to be fitted in the suite of ky models. Other numerical codes have been developed 
independently by several groups (Beckwith and Done [2004], C'adez and Calvani [2005], Brenneman and 
Reynolds [2006]) and equipped with similar functionality. 

However, the LAOR model can still be used for evaluation of the spin if one identifies the inner edge 
of the disc with the marginally stable orbit. In this case the spin is actually estimated from the lower 
boundary of the broad line. The comparison of the LAOR and kyrline model is shown in the right panel 
of Figure [T] The other parameters of the relativistic line models arc inclination angle i, rest energy of 
the line E, inner radius of the disc i?i n , outer radius of the disc i? ut, emissivity parameters qi, q± with 
the break radius r^. The emissivity of the line is given by / ~ r~ qi for r < and I ~ r~ q ' 2 for r > r^. 
The angular dependence of the emissivity is characterized by limb darkening profile J(/x e ) cx 1 + 2.06/i c 
in the LAOR model. The KYRLINE model enables to switch between different emission laws. We used 
further two extreme cases, the kyrline with the same limb-darkening law as in the LAOR model and 
KYRLINE* with the limb-brightening law I(fi e ) °c m(l + — ^). 

The aim of this paper is to compare the two models applied to the current data provided by the 
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XMM-Newton satellite, and to the artificial data generated for the on-coming X-ray mission. For this 
purpose we have chosen two sources, MCG-6-30-15 and GX 339-4, which exhibit an extremely skewed 
iron line according to recently published papers (Vaughan and Fabian [2004], Miller et al. [2004]). 
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Figure 1. Left: Relation for the spin a and marginally stable orbit r ms . Right: Comparison of the LAOR 
(black, solid) and kyrline (red, dashed) model for two values of the spin a = 0.9982 (top) and a = 0.7 
(bottom). The other parameters of the line are E — 6.4 keV, qi = qi = 3, i = 30°. 



Observations and data reduction 



We used the SAS software version 7.1.2 (http://xmm.esac.esa.int/sas) to reduce the XMM-Newton 



data of the sources. Further, we used standard tools for preparing and fitting the data available at 
|http://heasarc.gsfc.nasaiov| (FTOOLS, XSPEC) 

The galaxy MCG-6-30-15 is a nearby Seyfert 1 galaxy (z = 0.008). The skewed iron line has been 
revealed in the X-ray spectra by all recent satellites. The XMM-Newton observed MCG-6-30-15 for a 
long 350ks exposure time during summer 2001 (revolutions 301, 302, 303). The spectral results are 
described in Fabian et al. [2002]. We joined the three spectra into one using the ftool MATHPHA. 

The black hole binary GX 339-4 exhibited a strong broadened line in the 76 ks observation in 2002 
(Miller et al. [2004]) when the source was in the very high state (for a description of the different states 
see Remillard and McClintock [2006]). The observation was made in the burst mode due to a very high 
source flux. The 97% of photons are lost during the reading cycle in this mode, which results into 2.25 ks 
total exposure time0 

We rebinned all the data channels in order to oversample the instrumental energy resolution max- 
imally by a factor of 3 and to have at least 20 counts per bin. The first condition is much stronger 
with respect to the total flux of the sources - 4 x 10 -11 ergcm~ 2 s _1 in 2-10keV (1.1 x 10 6 cts) for 
MCG-6-30-15 and 9 x 10~ 9 ergcm" 2 s^ 1 in 2-10 keV (1.0 x 10 7 cts) for GX 339-4. 
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MCG-6-30-1 5: power law model r = 1 .9 MCG-6-30-1 5: Fe line 
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Figure 2. The X-ray spectrum of MCG-6-30-15 observed by XMM-Newton. Left: The overall view on 
the spectrum as a power law with r — 1.9 absorbed by a neutral hydrogen along the line of sight with 
column density tt-h = 0.41 x 10 21 cm -2 . The residuals from the model are plotted in the bottom panel 
clearly revealing features of a local absorption and soft excess at the soft X-ray band, and a feature at 
around 6kcV which can be explained by the presence of a broad iron line. Right: More detailed view of 
the iron line band. 



Iron line study of the MCG-6-30-15 spectrum 

We used the same continuum model for the MCG-6-30-15 spectrum as presented in Fabian et 
al. [2002]: the simple power law component absorbed by neutral hydrogen with the column density 
n.H = 0.41 x 10 21 cm~ 2 . The overall spectrum with a detailed view of the iron line energy band is shown 
in Figure^] The employed model is sufficient to fit the data above w 2.5 keV, which is satisfactory for 
our goal of the comparison of the relativistic line modelsd The value of the photon index is Y = 1.90(1). 
The residuals are formed by a complex of a broad iron line and two narrow iron lines - one emission line 
at E = 6.4 keV likely originating in a distant matter (torus) and one absorption line at E = 6.77 keV 
which can be explained by a blucshiftcd absorption originating in an outflow. The rest energy of the 
broad line is E = 6.7 keV, which corresponds to the helium-like ionized iron atomsH 

A good fit of the broad line was found with a broken power law line emissivity with a steeper 
dependence on the radius in the innermost region, which suggests a centrally localized corona. The 
goodness of the fit is constrained by the least squared method. The fit results in 2.5-9.5 keV are presented 
in Table 1. The y 2 values give comparable results for all employed models. The x 2 c d = X 2 l v ~ 1-2, 
where v is the number of degrees of freedom which is related to the total number of energy bins and 
model parameters. The six independent parameters of the LAOR and KYRLINE models make the global 
minimum of x 2 rather wide with several local minima. Each model has a different tendency to converge 
to a different minimum. Hence, we did not compare only best fits of both models, but also the evaluated 
spin values by the kyrline and LAOR models when the other model parameters correspond to each 
other. The equivalent width of the line is EW w 750 eV. The errors in brackets presented in the table 
correspond to 90% confidence and are evaluated when the other parameters of the model are fixed. The 
realistic errors are higher because the model parameters further depend on the other parameters of the 
line and continuum models. 

To catch up these relations we produce various contour graphs focusing on the determination of the 
spin value, taking into account the other parameters of the used model. The relations of the x 2 values 



1 The broad iron line was also revealed in the analysis of the two 138 ks observations in spring 2004 by Miller et al. 
[2006] , when the source was in the low-hard state. The EPIC pn camera was operating in the timing mode, MOS cameras in 
the full-frame mode. Using epatplot tool we found that it is not possible to avoid the pile-up by excluding the central part 
of the image of the source as described in the forementioned paper. Therefore we use only the very high state observation 
from 2002 in our analysis. 

2 Another components are needed to be added into the model in order to fully understand the spectrum. Several works 
have been done in this way, in the most recent one by Miller L. et al. [2008] the spectrum is characterized by an absorption 
in four different zones, which affects also the higher energy band where no broad line is needed any more. 

3 The spectral complexity in the line band allows an alternative explanation - the model with two narrow emission lines 
at energies E = 6.4 keV and E = 6.97 keV. This alternative model leads to the presence of the broad line component at 
E = 6.4keV. 
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on value of the spin (kyrline) or the inner disc radius (laor) are shown in the left column of Figure [31 
The x-axis is oppositely directed in the case of the inner disc radius as x- variable for an easier comparison 
with the kyrline results. The contour graphs for the spin and the inclination angle are shown in the 
middle column of Figure [H The underlying model was fixed in both cases. The plots in the right column 
of Figure [3j show the contours for the spin and the power law index. Taking all of these into account, we 
obtain for the spin value: 

a KY = 0.94I°;°6 and ^ = 0.96^8- 



Table 1. Results for MCG-6-30-15 in 2.5 9.5 keV 



parameter 


KYRLINE 


KYRLINE* 


LAOR best 


LAOR loc.min. 


a/M 
i [deg] 
E [keV] 

qi 

92 

n 


26.7(7) 
6.67(1) 
4.9(1) 
2.84(4) 
5.5(2) 


95+ - 02 
31.5(7) 
6.60(1) 
3.7(1) 
2.11(4) 
18.3(5) 


0.981°;°? 
35.7(5) 
6.48(1) 

4.8(1) 
2.50(3) 

6.6(2) 


96+° 02 

26.8(5) 
6.66(1) 

4.7(1) 
2.87(3) 

5.1(2) 


x l h 


175/148 


174/148 


170/148 


174/148 


EW [eV] 


761 


757 


764 


754 




spin a [GM/c] spin a/M [GM/c] spin a/M [GM/c] 



Figure 3. The contour graphs show the dependence of the value of the x 2 , the inclination angle and the 
power law index on the value of the spin (kyrline) or the inner disc radius (laor) for the MCG-6-30-15 
spectrum in 2.5-9.5 keV. The black, red and green contours correspond to ler, 2a and 3er, respectively. 
Top: The results of the LAOR model. Middle: The results of the kyrline model with limb darkening. 
Bottom: The results of the kyrline model with limb brightening. 
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Iron line study of the GX 339-4 spectrum 



GX 339-4: diskbb + powerlaw GX 339-4: Fe line 




Figure 4. The X-ray spectrum of GX 339-4. Left: The overall view on the spectrum as a power law 
with r = 3 and thermal multi temperature black body with 0.87 keV absorbed by a neutral hydrogen 
with column density nn = 0.6 x 10 22 cm -2 . The residuals from the model are plotted in the bottom 
panel revealing wiggles at the soft X-ray band which are likely due to instrumental response. The broad 
excess at around 6keV can be explained by the presence of a broad iron line. Right: More detailed view 
of the iron line band. 



The continuum of the X-ray spectrum of the black hole binaries is characterized by a power law 
and a multi-colour disc black-body component (powerlaw + DISKBB in the XSPEC notation). The 
power law index suggested from the simultaneous RXTE measurements is T = 2.5 (Miller et al. [2004]). 
However, we get an unacceptable fit with Xreducod — 6-7 using the same model as in the analysis by 
Miller et al. [2004] or the re-analysis by Reis et al. [2008] . The difference of the results is likely due to 
a different grouping of the instrumental energy channels applied to the data. While we did not allow to 
oversample the instrumental energy resolution more than by factor of 3, in the previous works only the 
condition to have at least 20 counts per bin was used. This condition is very weak with respect to the 
total number of counts N counts wl.Ox 10 7 and the total number of energy channels iV c han = 1.5 x 10 3 in 
2-10 keV and as a result, it practically does not force the data to be grouped. This leads to an excessive 
oversampling of the energy resolution, to large error bars in the flux and finally to an artificial decrease 
of the Xred vame - In the energy range 0.8-9 keV we get x 2 l v = 2835/1640 for the grouping with the 
only condition 20 counts per bin, and \ 2 l v = 1368/202 for the grouping taking the energy resolution 
into account. The x 2 e d value increased from x 2 e d = 1-73 to x 2 e d = 6.77. 

The absorption and emission features around E m 2 keV have a significant effect on the goodness of 
the fit. These features can be linked with the Si and Au edges suggesting to be instrumental imprints. 
Hence, we added two gaussian lines into the model to improve the fit (x 2 jv — > 667/202). As a next 
step, we allowed the continuum parameters to float (and also removed the SMEDGE model used in 
the previous works). The new model has x 2 l v = 350/202 in 0.8-9 keV and its parameter values are 
n H = 0.6 x 10 22 cm" 2 , T = 3.08(5) and kT in = 0.87(1) keV. We tried to add a reflection component 
into the model, as pexrav (pexriv) or refsch model, but without any improvements of the fit. The 
spectrum of GX 339-4 is shown in Figure 3] with a detailed view of the iron line band in the right panel. 
A broadened iron line feature is still present. However, due to different adopted value for the photon 
index of the power law the line is much weaker than the one presented in Miller et al. [2004] . 

The fitting results of the line models in 3-9 keV are summarized in Table 2 and Figure [5] There are 
two minima found during the fitting procedure. We preferred the one which better corresponds to the 
results obtained by the independent radio and infrared measurements which constrained the inclination 
angle to be i < 26° ( Gallo et al. [2004] ) . The dependence of the goodness of the fit on the spin value 
is shown in the left column of Figure [5] The contour graphs for the spin and the inclination angle are 
depicted in the middle column, and for the spin and the power law photon index in the right column of 
Figure [5] The derived spin value is then: 

a KY =0.691°;^ and a laor = 0.77^;!°. 
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Figure 5. The contour graphs show the dependence of the value of the y 2 , the inclination angle and the 
power law index on the value of the spin (kyrline) or the inner disc radius (laor) for the GX 339-4 
spectrum in 3-9 keV. The black, red and green contours correspond to la, 2a and 3tr, respectively. Up: 
The results of the LAOR model. Middle: The results of the kyrline model with limb darkening. Down: 
The results of the kyrline model with limb brightening. 



Table 1. Results for GX 339-4 in 3-9 keV Table 2. Count rates of the observations 



parameter 


KYRLINE 


KYRLINE* 


LAOR 


MCG-6-30-15 


GX 339-4 


a/M 
i [deg] 
E [keV] 

q 


0.89tS:S 

19(3) 
6.97(1) 
3.45(8) 


n R2+ ' 14 

U>D ^-0.14 

19(4) 
6.97(1) 
3.35(8) 


n --+0.08 
u> ' ' -0.12 

17(4) 
6.97(1) 
3.3(1) 


net cts/s 
model cts/s 
line cts/s 
line cts 


3.59 
3.59 
0.20 
4.37 x 10 4 


592.1 
592.5 
5.1 
1.15 x 10 4 


x7« 


147/125 


148/125 


148/125 








EW [eVJ 


175 


164 


164 









Fitting of the simulated data of future X-ray missions 

In this section we intend to apply the LAOR and kyrline models on the data with significantly higher 
quality supposed to be achieved by on-coming X-ray missions. The presently planned International X-ray 
Observatory (LXO) arised from the merging of the former XEUS and Constellation-X missions. Because 
the details of the IXO mission have not been fixed yet, we used a preliminary response matrix of the 
former XEUS mission (Arnaud et ai. [2008]). We generated the data for a kyrline model with a rest 
energy of the line E = 6.4 keV and radial disc luminosity that follows a power law with the index q = 3. 
We rebinned the data in order to have a resolution of 30 eV per bin. We then fit the data in the 1-9 keV 
energy range with the LAOR model using the same initial values of the fitting parameters as for the 
data simulation. Due to insufficient resolution of the LAOR model, a significant problem appears at the 
high-energy edge of the broad line. This occurs because the next generation instruments achieve much 
higher sensitivity in the corresponding energy range. Therefore, we excluded the higher-energy drop of 
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Figure 6. Simulated data for two values of the spin and inclination using the preliminary response 
matrix of the XEUS mission. We show the artificial data with the expected errors (black crosses) which 
were simulated by the KYRLINE model (black curve), and fitted by the LAOR model (red curve). 



the lines from the analysis in order to reveal the differences in the overall shape of the line. We examined 
the artificial data for a grid of values of the angular momentum and the inclination angle. The results 
for a/Al = 0.7, 0.9982 and i = 60°, 80° are shown in Figure[S] It is clearly seen from the figure that the 
effect of the spin on the shape of the line is sufficiently resolved by the higher quality data. 

Further, we produced simulated data for the Seyfert galaxy MCG-6-30-15 and the black hole binary 
GX 339-4 using rather simplified models which were suitable to fit the current XMM-Newton data. 
For MCG-6-30-15 we used a power law model plus a kyrline model for the broad iron line, absorbed 
by neutral hydrogen: phabs * (powerlaw + kyrline). The parameters of the continuum are the 
column density tt-h = 0.4 x 10 21 cm" 2 , the photon index T = 1.9 of the power law, and its normalization 
Ky = 5 x 10~ 3 . The values of the line parameters are summarized in the K Y value column of Table 4. 
The exposure time was chosen as 220 ks and the flux of the source as 1.5 x 10 _11 ergcm -2 s _1 in the 
2-10 keV energy range (i.e. 1.4 x 10 7 cts). 

For GX 339-4 we used phabs * (diskbb + powerlaw + kyrline) with n H = 0.6 x 10 22 cm~ 2 , 
fcT in = 0.87 keV (K kT = 1.4 x 10 3 ), and r = 3 (K r = 5.6). The exposure time was chosen as 75 ks, the 
flux of the source as 9.3 x 10 -9 ergcm~ 2 s _1 in the 2-10keV energy range (i.e. 1.5 x 10 9 cts)|f| 



4 The number of counts is two orders of magnitude higher than for the observation of the XMM-Newton satellite. The 
reason is due to the loss of 97% of the photons during the burst mode of XMM-Newton observation which eliminates 
the pile-up problem. The next generation X-ray missions are supposed to have a calorimeter instead of the CCD camera 
on-board which will get rid of such problems. 
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Prior to the spectral analysis we rebinned the data to have approximately a 5eV resolution (as it 
was planned for the XEUS instrument). We tested different grouping realizing that the discrepancies 
between the two models increase with larger grouping, but we can see apparent differences already for 
the most moderate rebinning. The results of the LAOR fit are shown in Table 4 and in Figure [7] The 
broad iron line component of the model is plotted in the left column. The continuum components of 
the model are not displayed there in order to clearly see the deflections of the LAOR model. The most 
prominent discrepancy appears at the higher-energy drop, which is clearly seen in the data/model ratio 
plot in the middle column of Figure [7] The model parameters are constrained with small error bars (see 
contours a vs. i in the right column of Figure [7|), which clearly reveals a difference between the kyrline 
and the LAOR models. 

The spin value derived from the analysis using the LAOR model is: 

— n QC.o+0.003 
fllaor.MCG — u - y O 8 _0.004 

for MCG-6-30-15, while initially the spin value was atrial, mcg = 0.940, and 

n 7/1+0.02 

aiaor,GX — U. 1 4_ Q Q2 

for GX 339-4, while initially the spin value was atrial, gx = 0.70. 



Table 4. Results of LAOR fit in 3 9 keV in the simulated spectra 





MCG-6-30-15 


GX 339-4 


parameter 


KY value 


fitted value 


KY value 


fitted value 


Rin [G] 
i[deg] 


2.00 
26.7 


1 Sf}+ - 04 
i - 8b -0.02 

25.4(3) 


3.39 
19 


o 20+ 05 

°- zu -0.05 

18.6(2) 


-Eline 


6.70 


6.71(1) 


6.97 


6.94(5) 




4.90 


4.51(3) 


3.45 


3.2(1) 


92 


2.80 


2.76(2) 






n 


5.5 


6.2(1) 






Kn ne 


8.7 x 10" 5 


9.0 x 10~ 5 


6.5 x 10~ 3 


6.6 x 10~ 3 


X 2 /v 




1364/873 




2298/873 



Conclusions 

We investigated the iron line band for two representative sources - MCG-6-30-15 (active galaxy) 
and GX 339-4 (X-ray binary). The iron line is statistically better constrained for the active galaxy 
MCG-6-30-15 due to a significantly longer exposure time of the available observations - for comparison 
of count rates of the sources see Table 3. The spectra of both sources are well described by a continuum 
model plus a broad iron line model. We compared modeling of the broad iron line by the two relativistic 
models, LAOR and kyrline. The kyrline model leads to a better defined minimum of % 2 for the best fit 
value. The confidence contour plots for a/M versus other model parameters are more regularly shaped. 
This indicates that the kyrline model has a smoother adjustment between the different points in the 
parameter space allowing for more reliable constraints on a/M. The LAOR model has a less accurate grid 
and is strictly limited to the extreme Kerr metric. The discrepancies between the kyrline and LAOR 
results are within the general uncertainties of the spin determination using the skewed line profile when 
applied to the current data. However, the results are apparently distinguishable for higher quality data, 
as those simulated for the XEUS mission. We find that the LAOR model tends to overestimate the spin 
value and furthermore, it has insufficient energy resolution which affects the correct determination of the 
high-energy edge of the broad line. The discrepancies in the overall shape of the line are more visible 
especially for lower values of the spin a/M . As a side-product, we have found that the correct re-binning 
of the data with respect to the instrumental energy resolution is crucial to obtain statistically the most 
relevant results. 
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Figure 7. The simulated spectra for MCG-6-30-15 (top) and GX 339-4 (bottom) using the preliminary 
response matrix of the XEUS mission. Left: The broad iron line generated by the kyrline model (black 
data) and fitted by the LAOR model (red curve). Middle: The data/model ratio when fitting by the LAOR 
model. Right: Contours for the inclination angle and the inner disc radius of the LAOR model. The trial 
values are far from the best fit results (see KY value in Table 4) . 



References 

Arnaud, M., Barcons, X., Barret, D., Bautz, M., Bellazzini, R., XEUS: the physics of the hot evolving universe, Experimental 
Astronomy, tmp, 24A, 2008 

Beckwith, K., Done, C, Iron line profiles in strong gravity, Monthly Notices of Royal Astronomical Society, 352, 353, 2004 
Brenneman, L. W., Reynolds, C. S., Constraining black hole spin via X-ray spectroscopy, Astrophysical Journal, 652, 1028, 
2006 

Cadez, A., Calvani, M., Relativistic emission lines from accretion discs around black holes, Monthly Notices of Royal 

Astronomical Society, 363, 177, 2005 
Dovfiiak, M., Karas, V., Yaqoob, T., An extended scheme for fitting X-ray data with accretion disc spectra in the strong 

gravity regime, Astrophysical Journal Supplements, 153, 205, 2004 
Fabian, A. C, Vaughan, S. Nandra, K., Iwasawa, K., Ballantyne, D. R. et al., A long hard look at MCG-6-30-15 with 

XMM-Newton, Monthly Notices of Royal Astronomical Society, 335, LI, 2002 
Gallo E., Corbel S., Fender R. P., A transient large-scale relativistic radio jet from GX 339-4, Monthly Notices of Royal 

Astronomical Society, 347L, 52G, 2004 
Laor, A., Line profiles from a disk around a rotating black hole, Astrophysical Journal, 376, 90, 1991 

Miller, J. M., Fabian, A. C, Reynolds, C. S., Nowak, M. A., Homan, J. et al., Evidence of Black Hole Spin in GX 339-4: 
XMM-Newton/EPIC-pn and RXTE Spectroscopy of the Very High State, Astrophysical Journal, 606, L131, 2004 

Miller, J. M., Homan, J., Steeghs, D., Rupen, M., Hunstead, R. W. et al., A Long, Hard Look at the Low/Hard State in 
Accreting Black Holes, Astrophysical Journal, 653, 525, 2006 

Miller, L., Turner, T. J., Reeves, J. N., An absorption origin for the X-ray spectral variability of MCG-6-30-15, Astronomy 
& Astrophysics, 483, 437, 2008 

Reis, R. C, Fabian, A. C, Ross, R., Miniutti, G., Miller, J. M., Reynolds, C. S., A systematic look at the Very High and 
Low/Hard state of GX 339-4: Constraining the black hole spin with a new reflection model, arxiv:0804-0238, 2008 

Remillard, R. A. and McClintock, J. E., X-ray Properties of Black Hole Binaries, Annual Review of Astronomy & Astro- 
physics, 49, 2006 

Vaughan, S., Fabian, A. C, A long hard look at MCG-6-30-15 with XMM-Newton - II. Detailed EPIC analysis and 
modelling, Monthly Notices of Royal Astronomical Society, 348, 1415, 2004 



